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We derive constraints on the asymmetry oi of the momentum fractions carried by quark 
and antiquark in K and K* mesons in leading twist. These constraints follow from exact 
operator identities and relate ai to SU(3) breaking quark- ant iquark-gluon matrix elements 
which we determine from QCD sum rules. Comparing our results to determinations of Oi 
from QCD sum rules based on correlation functions of quark currents, we find that, for 
a\{K*) the central values agree well and come with moderate errors, whereas for ai{K) and 
aj^{K*) the results from operator relations are consistent with those from quark current 
sum rules, but come with larger uncertainties. The consistency of results confirms that the 
QCD sum rule method is indeed suitable for the calculation of Oi. We conclude that the 
presently most accurate predictions for oi come from the direct determination from QCD 
sum rules based on correlation functions of quark currents and are given by: 

ai(X) = 0.06 ±0.03, a!(is:*) = 0.03 ±0.02, a^(/s:*) = 0.04 ± 0.03. 
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1 Introduction 



Hadronic light-cone distribution amplitudes (DAs) of leading twist have been attracting 
considerable interest in the context of B physics. They enter the amphtudes of QCD 
processes that can be described in coUinear factorisation, which include, to leading order 
in an expansion in 1/mf,, a large class of nonleptonic B decays [1], such as B ^ mijKK. 
DAs are also an essential ingredient in the calculation of weak decay form factors such 
as B ^ 71, p, K, K* from QCD sum rules on the light-cone [2]. These decays, and their 
CP asymmetries, are currently being studied at the B factories BaBar and Belle and are 
expected to yield essential information about the pattern of CP violation and potential 
sources of flavour violation beyond the SM. 

One particular problem in this context is the size of SU(3) breaking corrections to 
K and K* DAs, which has been studied in a number of recent papers [3, 4, 5, 6]. The 
DAs themselves are defined as matrix elements of quark-antiquark gauge-invariant nonlocal 
operators on the light-cone. To leading-twist accuracy, there are three such DAs for K and 
K* = 0): 



{0\q{z)^^,[z,-z]s{-z)\K{q)) = tMqz) / due'^^^'^K{u) , 

Jo 

{0\q{z)^[z,-z]s{-z)\K*{q,X)) = {e^'^z)flmK* f due'^^'^^Uliu), 

Jo 

\q{z)a,4z,-z]s{~z)\K*{q,X)) = i{efhu - e^^^ q,)fk{^) f due'^^^^m{u), (1) 

ig j daz'^Ai^[az) 



with the Wilson-line 

[z, — = Pexp 



inserted between quark fields to render the matrix elements gauge-invariant. In the above 
definitions, e^,^'' is the polarisation vector of a vector meson with polarisation A; there are 
two leading-twist DAs for vector mesons, and (pj^, corresponding to longitudinal and 
transverse polarisation, respectively. The integration variable u is the (longitudinal) meson 
momentum fraction carried by the quark, u=l—u the momentum fraction carried by the 
antiquark and ^ — u — u. The decay constants f^''^^ are defined in the usual way by the 
local hmit of Eqs. (1) and chosen in such a way that 

f du(p{u)^l. (2) 
Jo 

All three distributions (pK, ^i:, (pic can be expanded in Gegenbauer polynomials Cn'^, 

(t>{u) = 6«« [ 1 + ^ a^CTi'^u - 1) ) , (3) 

V n>l / 

where the a„ are hadronic parameters, the so-called Gegenbauer moments. 
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The most relevant quantities characterising SU(3) breaking of these DAs are the decay 
constants fx and /^'", and ai{K) and a^'^K*), which can be expressed in terms of the 
DAs as 

5 /"^ 

ai{K) = -j du{u-u)(j)Kiu) (4) 

and correspondingly for a\''^{K*). ai describes the difference of the average longitudinal 
momenta of the quark and antiquark in the two-particle Fock-state component of the 
meson, a quantity that vanishes for particles with equal-mass quarks (particles with definite 
G-parity). The decay constants fx and /|- can be extracted from experiment; has been 
calculated from both lattice [7] and QCD sum rules, e.g. Ref. [6]. In this paper we focus on 
the determination of ai. no lattice calculation of this quantity has been attempted yet, so 
essentially all available information on ai comes from QCD sum rule calculations, ai can be 
calculated either directly from the correlation function of two quark currents [3, 4, 6, 8, 9] or 
from operator identities relating it to certain quark-quark-gluon matrix elements, denoted 
^4, which are calculated from QCD sum rules themselves [5]. In a previous paper, Ref. [6], 
we have obtained the following results from the first method, at the scale of 1 GeV: 

ai(/s:)^2 = 0.050 ±0.025, a!(/s:*)^^ = 0.025 ± 0.015, a^(ir*)^^ = 0.04 ± 0.03, (5) 

whereas Braun and Lenz found the following results from operator identities [5] : 

ai{Kf^ = 0.10 ±0.12, al(ii'*)^^ = 0.10 ±0.07. (6) 

These results were obtained to first order in and neglecting explicit terms in and 
rriq in the operator identities. Numerically, however, these terms are not negligible: the 
O(m^) correction shifts ai{K) by +0.17 and a\{K*) by +0.08 for our central value of m^. 
Corrections in rriq are small for a\{K*), but chirally enhanced for ai{K) and shift ai{K) 
by +0.04 for our central value of niq. A consistent inclusion of 0{mq^s) effects requires the 
calculation of these terms also for K4. In the present paper, we present such a calculation 
and improve the sum rules for K4 derived in Ref. [5] by the inclusion of all dominant terms to 
0{mq) and 0{ml), which include in particular two-loop perturbative and gluon-condensate 
contributions. The perturbative contributions come with large coefficients and prove to be 
very relevant numerically. We then construct several sum rules for K4 which differ by the 
chirality structure of the involved currents and the spin-parity assignment of the hadronic 
states coupling to them. We provide criteria that allow one to identify the sum rules most 
suitable for the calculation of and obtain the corresponding numerical results, including 
a careful analysis of the theoretical uncertainty of K4 and the corresponding values of ai. 
One important finding of our paper is that the results of these calculations agree, within 
errors, with those from the quark current sum rules, which shows that the application of 
the QCD sum rule method to the calcualation of oi yields mutual consistent results. It is 
this consistency that strengthens our confidence in the validity of the results for oi. 

Our paper is organised as follows: in Sec. 2 we derive the operator relations for oi, 
in Sec. 3 we obtain numerical results for the corresponding matrix elements and compare 
with the results of Ref. [6] . In Sec. 4 we summarise and conclude. The paper also contains 
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two appendices giving explicit expressions for all relevant correlation functions and Borel 
transforms. 



2 Exact Identities for ai 

In Ref. [5], the following relations were obtained for ai{K) and a\{K*): 

5 = + ' - '"^^"^^ ' 

\ al{K*) = -4 + 2 - 44(i^*), (8) 

5 f\\^ rriK* mj^, 

where K4^{K) and k\{K*) are twist-4 quark-quark-gluon matrix elements defined by 

{0\q(gG^i,)iri5s\K(q)) = iqjKmlK^(K), (9) 

{0\q{gG^,)trs\K*{q, A)) = eW/im^,4(i^*). (10) 

K4(K) and K4(i^*) vanish for — > due to G-parity. The special structure of (7) allows 
one to determine the value of K4^{K) to leading order in for niq [5], 

>^4iK) = (11) 

which is a consequence of the conservation of the axial current in the chiral limit. 

The above relations were derived from the analysis of matrix elements of the local 
operators {D—D — D) 

1 ^1 ^1 ^ 

of} = 2 ^lni.l^)i DvS + - qiu{lb)i D^^s-- g^,^qi{^5) P s, (12) 

whose divergence can be expressed in terms of bilinear quark operators. In this section, 
we rederivc these relations in a different way and obtain a new relation for ai{K*). 
The starting point for our analysis are the exact nonlocal operator relations [10, 11] 



d 



/•I 

- q{x)1h{15)s{-x) = -I j dvvq{x)xagG'^''{vx)j^{-f5)s{-x) 

M J-1 

- {ms ± mq)q{x)i{^^)s{-x), (13) 



d^{(l{x)^n{lb)s{-x)} = - ij dv q{x)xaG"''{vx)'y^{'y5) 



s[-x 



{mqTms)q{x)i{'j5)s{-x), (14) 



where the total translation is defined as 



d 

{q{x)Ts{-x)] = — {q{x + y)[x + y, -x + y]Ts{-x + y)} 



(15) 



The corresponding nonlocal matrix elements are, for K and K*^ {x"^ ^ 0): 

{0\q{xh^,j,s{-x)\K{q)) = i^q^ f due'^^^ + 0{x')] 

Jo 

+ ^ Um^K — du e^^«^ [gK{u) - M^) + ^(x^)] , (16) 

^ qx Jo 

mx)t^,s{-x)\K{q)) = -^^^ fdue^^'^^ + , (17) 

777.5 + T^q Jo 

II { e^^^x r II 1 

m{x)^,s{-x)\K*{q,X)) = flruKA q, / due'^'^' <t>l{u) + 0{x^)\ 

\ qx Jo '- 

+ - g^^) due^^'^- {gl{u) + 0(^2)) 

v qx / Jo 



In the above definitions, (pK and 0^ are the leading- twist DAs of K and Kp respectively; 
all other functions are higher-twist DAs and have been extensively discussed in Refs. [10, 
11, 12, 13]. 

ai{K), the quantity we are interested in, is related to the first moment of (pK' 

a,{K) = ^-Mt 

with M/ = J(^^ du {u — u)f{u) being the first moment of the DA f{v). Taking the matrix 
elements of (13) and (14) for K and expanding to leading order in x^ and next-to- leading 
order in qx, one obtains the exact relations 

rus + niq 

^ (^Mf^ + Mf^) = -2K4iK) + Mf", (19) 

from which one can determine M{ ^ once either ^ or Mf ^ are known, qk is a twist-4 
DA and Mf^ contains quark-quark-gluon matrix elements itself, cf. Refs. [11, 13], whereas 
0^ is twist- 3 and M-^'^ is completely determined in terms of the twist- 2 DA (t)K and 
mass corrections. ^ can be obtained from a second set of nonlocal operator relations 
involving tensor currents q{x)a ^jj'^^s{—x) or, equivalently, from the recursion relations for 
the moments of 0^ given in Ref. [13]: 

xp mj. — ml 
Ml" = , ^ 
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Solving (19) for ai{K), we then rederive 
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nis + niq 



1 1 
mz — m,. 



m 



(20) 



K 



which confirms the result obtained in Ref . [5] . Note that the first term on the right-hand 
side is rather sensitive to the value of niq and the second one to that of m^. 
For Kp the same method yields the equations 



nis - m, 



Ml'' - Ml 



(3) 

9k' 



1 2 

mi — mi II 



m 



(21) 



K* 



Again, Qj^ is a twist-4 DA whose first moment is not known from any independent analysis, 
whereas Mf^, the first moment of the twist-3 DA g]^, can be read off Eq. (4.6) in Ref. [12]: 

(22) 

Ik 

We can then solve (21) for a\{K*) and obtain 

2 2 

m — m 

(23) 



5 jll, mK* 



which agrees with Eq. (8), the result obtained in Ref. [5]. 

Let us now apply the same method to chiral-odd operators, with the aim of obtaining 
an analogous new expression for ai{K*). The relevant nonlocal operator relations are 

d 



dx, 



-q{x)a^„s{-x) = -idt,q{x)s{-x) + (m^ - mq)q{x)'-f„s{-x) 



+ J dvq{x)gXaG"j^{vx)s{—x)—i J dvvq{x)gXaG'^'^{vx)afj,^s{—x), 



d ^ 

d^{q{x)a^i,s{-x)} = -i-— q{x)s{-x) - {ms + mg)q{x)-f„s{-x) 

ox,, 

+ / dvvq{x)gXaG"^{vx)s{—x)—i / dvq[x)gXaG'^^{vx)a^ys[—x). 
7-1 J-i (24) 

These relations were first derived, without the terms in ±mq, in Ref. [10]; the terms in 
the quark masses are new. 

The relevant K* matrix elements are given by [10]: 



\q{x)a,,s{-x)\K*{q,\))=zfk 



Jo 



<f>j,{u) + 0{x') 



5 



(qx) 



(25) 



|g(a;)s(-a;)|i^*(g,A)) 



K 



rriK* 



m 



K* 



due'^'i'^ + 0(a;2)) ,(26) 



where, again, 0^ is the leading-twist DA of the transversely polarised K* and /i^* and h^^ 
are higher- twist DAs. In addition, we also need the following quark-quark-gluon matrix 
element: 

m{9GJ')ap,s\K*{q,\))^ 

- + e^Q.) + ^i{K*){e^^\p - e^^^gjj . (27) 

Here ^3(7^'*) is a twist-3 matrix element, k-^{K*) is twist-4; both are 0{ms — ruq) due to 
G-parity.^ Taking matrix elements of (24), one obtains expressions in q^. e"u^ and x^,. To 
twist-4 accuracy only the former two are relevant and yield a set of four linear equations 
for the four first moments of gf^, /i^, /i^ and h^^: 



kUk*) + 2Ki{K*) + S+Mf"" - Mf^ - 



,(3) 



2^2^' 



1 



,(3) 



1 



l(3) 



3 m;^ - ^ = 25_ , 
- 



(28) 



with 5± = 7T ■ The solution of that system implies 



0" 1 
^ ^ 6 



i(5+Mf- + lMf--24(i^*), 



which must agree with Mf^ as given in Eq. (22). Solving for ai{K*), one finds 

which is the wanted new relation for a^{K*). Note that in all three relations (7), (8) and 
(29) enters multiplied with a large numerical factor which implies that the theoretical 
uncertainty of the resulting values of ai will be much larger than that of itself. 

^The normalisation of k^{K*) is chosen in such a way that / 'DQj'{q) = k^{K*) for the twist-3 DA 
T(a) defined in Ref. [12]. 
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{qagUq) = 


(-0.24 ± 0.01)^ GeV^ 

9 / — \ 

^0 m) 


(ss) = (1 - Ss) (qq) 
{sagGs) = (1 - d5){qagGq) 




(0.012 ± 0.003) GeV^ 




ml = 


= (0.8 ± 0.1) GeV^ 5^ 


= 0.2 ±0.2, = 0.2 ±0.2 


m,(2GeV) 


= (100 ± 20) MeV < — > m,(l GeV) = (137 ± 27) MeV 

rriqiiJ,) =ms{iJ,)/R, i? = 24.4 ±1.5 


as{mz) - 0.1187 ±0.002 ^ 


a,(lGcV) - 0.534+°-°^^ 


fK = 


(0.160 ± 0.002) GeV, /| = (0.217 ± 0.005) GeV 
= (0.185 ± 0.010) GeV 



Table 1: Input parameters for sum rules at the renormalisation scale /i = 1 GeV. The value 
of nis is obtained from unquenched lattice calculations with n/ = 2 flavours as summarised 
in [14], which agrees with the results from QCD sum rule calculations [15]. friq is taken 
from chiral perturbation theory [16].^ Ois(mz) is the PDG average [18], fx and /|- are 
known from experiment and has been determined in Refs. [6, 7]. The errors of quark 
masses and condensates are treated as correlated, see text. 

3 QCD Sum Rules for A and hij- 

In order to obtain numerical predictions for ai from the relations derived in the last section, 
one needs to know the values of the ^4 matrix elements. K4,{K) and kI{K*) have been 
calculated in Ref. [5] from QCD sum rules to leading order in SU(3) breaking parameters 
with the following results: 

K4(ir)^L = -0.11 ±0.03, = -0-050 ±0.010, (30) 

which, using the relations (7) and (8), letting = and neglecting the terms in 
translates into [5] 

ai(X)^L = 0.10 ±0.12, a!(X*)^^ = 0.10 ±0.07. (31) 

All these results refer to a renormalisation scale of 1 GeV. 

In this section we present QCD sum rules for ^^(K) and kI(K*) which are accurate to 
NLO in SU(3) breaking and also a new sum rule for (X*) to the same accuracy. For all 
sum rules we include 0{m^ effects. The sum rules are of the generic form 

«4(ir)/|-m^e"""*^'^^^ ± contribution from higher mass states = Bm'^^g-, (32) 

^rug has also been determined from lattice calculations. The most recent papers on this topic are 
Refs. [17]. The central value of mg/mq determined in the first of these papers with n/ = 2 running flavours 
and nonperturbative renormalisation agrees with the result from chiral perturbation theory, whereas the 
result of the second, obtained with n/ = 3 and perturbative (two- loop) renormalisation, is a bit lower. As 
the field appears to develop rapidly, we refrain from taking either side and stay with the result from chiral 
perturbation theory. 



7 



and correspondingly for K* . Hq are correlation functions of type 

with suitably chosen Dirac structures and explicit expressions for all relevant 
are given in App. A. Bm'^^g is the Borel transform of Uq, M'^ the Borel parameter and 
n is either 2 or 4. In order to separate the ground state from higher mass contributions, 
one usually models the latter, using global quark hadron duality, by an integral over the 
perturbative spectral density: 

contribution from higher mass states ~ / e~^^^ — ImnG(s); (33) 

the parameter sq is called continuum threshold. The input parameters for the QCD sum 
rules are collected in Tab. 1. 

All parameters can actually be determined from more than one sum rule derived 
from various which can be characterised by the following features: 

• the currents can have the same or different chirality, which results in chiral-even and 
chiral-odd sum rules, respectively; 

• the hadronic states saturating Hq can have unique spin-parity or come with different 
parity (e.g. 0~ and 1"*"), which results in pure-parity and mixed-parity sum rules, 
respectively. 

Note that all chiral-odd sum rules are also pure-parity. 

In chiral-odd sum rules the quark condensates always appear in the combination (gg) — 
{ss) = 5^{qq) and (qcrgGq) — (sagGs) = 5^{qagGq) , which induces a large dependence 
on the only poorly constrained parameters ^3^5 and also increases the impact of the gluon 
condensate contribution which is equally poorly known. We therefore decide to drop all 
chiral-odd sum rules and only use chiral-even ones. 

As for mixed and pure-parity sum rules, they come with different mass dimensions: 
n = 2 in (32) for mixed-parity vs. n = 4 for pure-parity sum rules. It is an important 
result of this paper that the continuum contributions to the mixed-parity sum rules, for 
typical Borel parameters around 1.7 GeV^, are small and below 10% for all three k^. 
Pure-parity sum rules, on the other hand, have a large continuum contribution around 30%. 
There are two reasons for this result: first, the additional power of in pure-parity sum 
rules counteracts the exponential suppression of the continuum contribution. Second, the 
contributions of particles with different parity have different sign: it was already found in 
Ref. [5] that H/i{K) and k}\{Ki) have opposite sign; we find that the same applies to ^{(yK*) 
and the corresponding n^{Kl) of the lowest scalar resonance, and ditto to k^{K*) and the 
couphng K^{Ki) of the axial vector Ki meson. These results suggest that the matrix 
elements of opposite-parity mesons have generically different signs and tend to cancel each 
other in mixed-parity sum rules, which results in a small continuum contribution. From 
a more formal point of view it is rather obvious from the definitions Eqs. (9), (10) and 
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(27) that the sign of K4 changes under a parity transformation,^ which is in hne with our 
findings. 

The mixed-parity sum rules for K and K* do involve the three spin-parity systems 
(0~, 1+), 0"*") and (1", 1+). Note that for all of them the "wrong" -parity ground state 
(e.g. the scalar Kq{1A30)) and the first orbital excitation of the "right" -parity state (e.g. 
the vector ir*(1410)) have nearly equal mass, which makes the cancellation very effective. 
We conclude that mixed-parity sum rules are more reliable than pure-parity ones and, as 
a consequence, will not consider the latter in this paper. In view of the cancellation of 
contributions of different sign we also decide to include explicitly only the lowest-mass 
ground state in the mixed-parity sum rules, which differs from the procedure adopted by 
the authors of Ref. [5]. 

Let us now turn to the question how to choose the Borel parameter and the contin- 
uum threshold sq, the internal sum rule parameters. As mentioned before, the dependence 
of the sum rules on sq is weak and so we simply use the same values of sq as for the 
quark current sum rules, i.e. So{K) = (1.1 ± 0.3) GeV^ 4{K*) = (1.7 ± 0.3) GeV^ and 
s^{K*) = (1.3 ± 0.3) GeV^ [6]. The small dependence on sq also allows one to use slightly 
higher values of than the usual 1 to 2GeV^, which improves the convergence of the 
operator product expansion of the correlation functions and reduces the variation of the 
sum rule with M^. We choose = (1.6 ± 0.4) GeV^ for /sT and = (1.8 ± 0.4) GeV^ 
for K*. 

After this general discussion of the choice of sum rules and parameters let us now turn 
to the three K4 parameters in turn. 

3.1 H,4{K) 

The mixed-parity sum rule for K4[K) is obtained from the correlation function 11^2 in 
App. A, Eq. (A. 6), and given by 

flmU^K) e-^/^^ = ^ - [ ds e^/^^ (lO In - 25 

+ Q — {rn,{qq) - m^{ss)) <^ -- + 7^ - In — + / — e"^/*^ 

1 Oi 1 

+ -g- ^ (^s(ss) - mq{qq)) + {ms{sagGs) - mq{qagGq)) 

+ ^ 1(99)^ - ■ (34) 

This sum rule includes all relevant contribiitions up to dimension six. Numerically, all 
dominant contributions have the same sign, with the largest one from (ss) , followed by the 
^In QCD parity is not a symmetry of the hadronic spectrum because the [/'(l)A-symmetry is broken. 
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-0.06 - K4{K) 




1-2 1.4 1.6 1.8 2.^/2 1 _ 2 1.4 1.6 1.8 2 . 

Figure 1: (Colour online) Left panel: K4^{K) from (34) as function of the Borel parameter 
M^. Parameters: renormalisation scale /i = 1 GeV, sq = 1.1 GeV^. Solid line: central value 
of input parameters; dashed lines: red: ^3 = 0,0.4, green: ^5 = 0,0.4, blue: ^^(m^) = 
0.1167, 0.1207. Right panel: ai{K) as function of from the operator relation (7) (colour- 
coded as in the left panel) and the value of ai{K) determined in Ref. [6] (purple lines). 
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Figure 2: Dependence of the central value of a^iK^ from (7) on 771^(2 GeV). Solid line: 
7715(2 GeV) = 4MeV, dashed line: 7775(2 GeV) = OMeV. 



ones from (sagGs) and perturbation theory which are roughly of the same size. 

In Fig. 1 we plot the resulting values for K4{K) and, via (7), ai{K), displaying, for 
illustration, explicitly the dependence on as and 53^5. It is evident that the dependence of 
both quantities on ^3 and ^5 is nonnegligible; at the same time, the comparison with ai{K) 
obtained in Ref. [6] from a QCD sum rule for quark currents shows that both sum rules 
agree within errors. Note that the inclusion of the perturbative contribution is crucial: 
without it, we would have obtained a negative result for ai{K). The impact of nonzero 7^7^ 
is also relevant and shifts the central value of ai{K) by +0.025. 

As for the theoretical uncertainties of K4{K) and ai{K) we note that they arise first from 
the QCD sum rule parameters and second from the uncertainties of the hadronic parameters 
given in Tab. 1. As for the former, as already stated above, we choose = (1.6±0.4) GeV^ 

*The results from the quark current sum rules quoted in this paper are slightly larger than the ones given 
in Ref. [6]. This is due to the fact that we have included infrared sensitive terms of type m'^ \n{M'^ / rrig) 
in the contribution of the gluon condensate in the mixed quark-quark-gluon condensate rather than in the 
Wilson-coefficient of the gluon condensate, cf. the discussion in App. A and Ref. [19]. 
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and So = (1.1 ± 0.3) GeV^ and add the corresponding uncertainties in quadrature. As for 
the latter, we treat rrig, rris, (qq), {(i<^gGq), 5^ and ^5 as parameters with correlated errors. 
Chiral perturbation theory helps to unravel some of these correlations: for instance, one 
has {ms + mq) / {2mq) = mj^/ml and m|- = —2{ms + mq){qq)/fl in LO chiral perturbation 
theory [16]. The dependence of 53,5 on is unfortunately unknown (and indeed would 
deserve further study). In order to estimate the uncertainty of Ki{K) and ai{K), we hence 
eliminate, using the above relations, rUq and {qq) as independent parameters in favour 
of iris-, but keep mg = (qcrgGq) / (qq) and ^3^5. This procedure is likely to overestimate 
the uncertainties induced by (ss) and (sagGs), but it is difficult to do better at present. 
Varying all remaining independent input parameters within their respective ranges given 
in Tab. 1, we obtain the following results: 

K^(K) = -0.09 ± 0.01 ± 0.01 ± 0.01 ± 0.02 ± 0.01 ± 0.00 = -0.09 ± 0.01 ± 0.02, 

ai{K) = 0.07 ± 0.04 ± 0.03 ± 0.11 ± 0.07 ± 0.03 ± 0.01 = 0.07 ± 0.04 ± 0.14, (35) 

where the first uncertainty comes from the variation of the sum rule specific parameters 
and sq, the second one from ag, the 3rd from m^, the 4th from ^3, the 5th from ^5 and 
the 6th from itLq = (qcrgGq) / (qq) . For the total uncertainty we give two terms: the ffist 
comes from the sum rule parameters and the second is obtained by adding all hadronic 
uncertainties in quadrature. As mentioned before, any uncertainty of k.4{K) induces a 
corresponding uncertainty in ai{K) that is about four times larger, except for the strange 
quark masses whose uncertainty also plays in the second term on the right-hand side of (7). 
The dependence of ai{K) on rus is shown in Fig. 2. The effect of nonzero rUg in the first 
term on the right-hand side of (7) is a shift by -1-0.04, which is partially, but not completely, 
compensated by the mg-dependent contributions to n^i^K). Comparing with the value of 
ai{K) quoted in Ref. [5], Eq. (6), we see that the central value in (35) is smaller and also 
the total uncertainty is larger. The larger error is due to the fact that we have chosen 
slightly larger errors for and also have included the uncertainty induced by a^. 

Let us now compare the result (35) with the one obtained from quark current sum rules 
[6], with the same sequence of errors as in (35): 

ai{Kf^ = 0.06 ± 0.01 ± 0.00 ± 0.01 ± 0.01 ± 0.01 ± 0.00 = 0.06 ± 0.01 ± 0.02 . (36) 

This number is slightly larger than the one quoted in Ref. [6], cf. footnote 2. Although the 
central values of ai{K) agree very well and hence confirm the consistency of the sum rule 
results, it is obvious that the operator relation (7) cannot match the accuracy of the quark 
current sum rule and is hence not very useful for constraining ai{K). 

3.2 kI{K*) 

Let us now turn to k'|(/C*). The mixed-parity sum rule is derived from the correlation 
function IIq\ ^PP- ^1 Eq. (A. 11), and reads 

4miArmU-<^^^' = {ml - ml) ^ £ dse-^'^' (lOln - 25^ 
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Figure 3: Left panel: i^\{K*) from (37) as function of the Borel parameter M^. Parame- 
ters: renormalisation scale jj, = 1 GeV, sq = 1.7 GeV^. Solid black line: central values of 
parameters; the coloured lines have the same meaning as in Fig. 1. Right panel: a\{K*) 
as function of from the operator relation (8) and the sum rule for a\{K*) calculated in 
Ref. [6] (purple lines). 



2«., . , , /- \^ / 1, . M\ f^ds 

(m,(gg) - mg{ss)) <^ -- + 7^5 - In — + / — e 



13 fJ' Jso S 

+|^i{OT>^-t«n. (37) 

The resulting values of k\{K*) and a\{K*) are shown in Fig. 3. Again, the contribution 
from perturbation theory is crucial numerically: without it, the resulting values of a\{K*) 
would have been negative. Our final results are: 

kI{K*) = -0.022 ± 0.003 ± 0.001 ± 0.003 ± 0.004 ± 0.001 ± 0.001 



= -0.022 ±0.003 ±0.005, 

al{K*) = 0.01 ± 0.02 ± 0.01 ± 0.01 ± 0.02 ± 0.01 ± 0.00 

= 0.01 ±0.02 ±0.03 (38) 

with the same assignment and treatment of uncertainties as in (35); the uncertainty coming 
from is included in that from mg. In contrast to the pseudoscalar case, the translation 
of k\{K*) into a\{K*) does not increase the uncertainty from any more than the other 
uncertainties, so that the total error of a\{K*) is smaller than that of ai{K). The impact 
of m^-dependent terms in negligible. The results (38) differ from those of Ref. [5], (30) and 
(31), where the pure-parity sum rule has been used instead. The result from the quark 
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Figure 4: Left panel: Kj;{K*) from (40) as function of the Borel parameter M^. Parame- 
ters: renormalisation scale /x = 1 GeV, Sq = 1.3 GeV^. Solid black line: central values of 
parameters; the coloured lines have the same meaning as in Fig. 1. Right panel: ai{K*) 
as function of from the operator relation (29) and the sum rule for aj^{K*) calculated 
in Ref. [6] (purple lines). 



current sum rule is 



al{K 



*\BZ 



0.03 ±0.02. 



(39) 



Again we find agreement between the results for ai from the sum rules for ^4 and the quark 
current sum rules, but at the same time the uncertainty of the former is larger than that 
of the latter. 



3.3 

The last parameter left to be determined is Kj;{K*). Its mixed-parity sum rule is derived 
from the correlation function nG,4, Eq. (A. 16), and reads 



a. 



727r3 



rfse"^/^^' f -61n-^ + 14 



+ 



^liQQ) - 2(ss) j - ^^?^|^ {ss) - 2(gg)| + {mq{qagGq) - m,{sagGs)) 



3n 



37r 



(40) 



The results for k,^{K*) and aj^{K*) are shown in Fig. 4; including uncertainties, we find 
nj{K*) = 0.018 ± 0.004 ± 0.001 ± 0.002 ± 0.002 ± 0.002 ± 0.001 

= 0.018 ±0.004 ±0.004, 
ai:{K*) = 0.09 ± 0.04 ± 0.01 ± 0.01 ± 0.02 ± 0.02 ± 0.01 

= 0.09 ±0.04 ±0.03. (41) 
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Note that the "enhancement" factor of uncertainties of ai{K*) due to /tf (fC*) is 10, which 
is the reason for the large total uncertainty in (41). The impact of m^-dependent terms is 
again negligible. The quark current sum rule yields [6] 

ai{K*)^^ = 0.04 ± 0.01 ± 0.01 ± 0.01 ± 0.01 ± 0.00 ± 0.00 = 0.04 ± 0.01 ± 0.02. (42) 

Hence, also for aj^{K*) do the results of the two approaches agree within errors, with the 
quark current sum rule being more accurate. 



4 Summary and Conclusions 

In this paper, we have obtained the following relations for the first Gegenbauer moments 
of the leading-twist distribution amphtudes of K and K* mesons: 

9 , . ms-rriq , ml-m\ 

- ax\K) = ■ ^ + 4 2 — - ^i^4{K), 

5 rUs + rUq mj^ 

5 ' ' /I rriK* m|f* 



where the K4 matrix elements are defined as 

{0\q{gGa^^)iY15s\K{q)) = iqafK'm\Ki{K) , 
{m9G.^Ws\K*{q,\)) = eW/l m|,.4(X*) , 
m{9GX)<y0,s\K\q,X)) = 

= \]^4{K*){e^ahp + efq^) + 4(^*)(eL'^g/. - efq^^ . 

The first two relations in (43) were already derived in Ref. [5], the third is new. We have 
interpreted these relations as constraints on oi and calculated the three parameters from 
QCD sum rules. We have improved the sum rules given in Ref. [5] for ki{K) and k\{^K*) by 
including two-loop pcrturbativc contributions, the gluon condensate contribution and terms 
in mq\ the former proved to be very relevant numerically, the terms in are relevant for 
ai{K). We have also derived a new sum rule for k-^{K*) to the same accuracy. All these sum 
rules exhibit only a small continuum contribution and all relevant contributions come with 
equal sign. The results for oi obtained from the relations (43) agree, within errors, with 
those obtained in Ref. [6] from quark current sum rules which is an important confirmation 
of the consistency of QCD sum rule calculations of these quantities and strengthens our 
confidence in the results. Prom a phenomenological point of view, however, the operator 
relations (43) are, at least at present, less useful than the quark current sum rules for oi, as 
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the uncertainties of the K4 parameters are too large to allow an accurate determination of 
Oi- The uncertainties of K4 arise from (a) the dependence of the sum rule on the sum rule 
internal parameters and sq, (b) the uncertainties of at the hadronic scale ~ 1 GeV 
and (c) the uncertainties of nis and the SU(3) breaking of quark and mixed condensates 
parametrised by 63^5. All these uncertainties enter oi multiplied by large factors 5 to 10, 
Eqs. (43). In contrast, the quark current sum rules for ai studied in Refs. [4, 6] are not 
very sensitive to these effects and come with smaller uncertainties. We hence suggest that 
the relations (43) be interpreted as constraints on rather than ai. Using the updated 
values of oi from quark current sum rules quoted in Sec. 3, adding the errors linearly, 

ai(X)B2 = 0.06 ±0.03, a! (X*)'^^ = 0.03 ±0.02 a^(X*)'^^ = 0.04 ± 0.03, (44) 

we find by solving (43) for K4: 

, , 1 iris - rrio 9 , , ml - ml 
kAK) = -- — ^ - — aAK) + I ^ = -0.09 ± 0.02, 

ms-mg 3 ,|,^^,, , m^-ml 



kI(K*) = ' ' - — a\{K*) + ^ , ^ = -0.024 ± 0.003, 

f l ^rriK* 20 ' 2mi. 



For K4{K) and k\{K*) the central value agrees well with the results from the direct calcu- 
lation, for Kj-{K*) there is agreement within errors. How can these results be improved? 
The quark current results for ai would profit from a calculation of perturbative radiative 
corrections ~ m^cts, which is technically feasible, but beyond the scope of this paper. Both 
ai and K4 would benefit from a reduction of the errors of m^. 

In summary, we hope that the present paper helps to settle the controversy about ai 
which started from the observation that the original calculation of Ref. [8] suffers from a 
sign-mistake of the perturbative contribution, which was corrected in Ref. [3]. Unfortu- 
nately, the chiral-odd sum rules used in Ref. [3] come with large cancellations of the domi- 
nant contributions and are hence not very useful for precise calculations of ai. In Ref. [4], 
ai{K) was then determined from chiral-even quark current sum rules and in Ref. [6] also 
a[^'^^\K*) was calculated using that method. These sum rules do not exhibit any cancel- 
lations of large contributions and are stable under the variation of all input parameters. 
As we have shown in this paper, these results agree with those from the operator relations 
(43) within errors, but are more accurate. We conclude that the quark current sum rule 
results (44) present the presently best determination of Oi. Given the phenomenological 
importance of ai, an independent calculation on the lattice would be both timely and useful 
and we would like to appeal to the lattice community to take up the challenge. 
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A Correlation Functions 



In this appendix we give the relevant formulas for the correlation functions from which the 
QCD sum rules given in Sec. 3 are obtained. The correlation functions are of the generic 
form 

n„...(g) = i I d'ye'^y{0\T[q{gG^,)r^,s]{y)[sr2q]{0m , (A-l) 

where and r2 are suitably chosen Dirac structures. The dots stand for additional indices 
from 



A.l K4{K) 

K.4{K) can be extracted from either a pure-partity sum rule, to which only pseudoscalar 
states contribute, or a mixed-parity sum rule which also contains contributions from axi- 
alvector mesons. As for pure-parity sum rules, one possible choice of the Dirac structures 
is = and r2 = ^75, which results in the correlation function 

n„(g) =^g„ng)(g2). (A.2) 

Another choice is = i'y'^'y^ as before and r2 = 7/375, with the correlation function 

liafsiq) = gaf^^l^\{q') + qaqpIi^oU^'): (A.3) 

where receives contributions from 1+ intermediate states only, whereas 11^2 is a 
mixed-parity correlation function with contributions from both 0~ and l"*" states. 

These three correlation functions are not independent of each other, but related by 

d^s^plbq — i'lTT's + ^q)si'y5q, SO that 

n^\(g2) + g^n^^^2(?^) = i^s + mq)u'^^ + contact terms, (A.4) 

where the contact terms are independent of q^. As terms in niq are numerically relevant in 
the operator relation (7), we calculate the correlation functions to 0{mq) and find 



iq') 



a. 



- m., - m„ 



487r3 



In^ 



In- 



Aq'^ 



[{qagGq) - {sagGs)] 



+ 



^[{qq)-{ss)]\n- ^ 
1 



8g2 



a. 



TT 



G' 



1 - In 



— I 1 — In 



/ 2 2\ 



s 2 

q 



71n^ 



-q 



11^ 



-471n- 



Stt 



-{qq) - (ss) 



In- 



-q^ ruqas 
37r 



{qq) 



In 



rriq ms \ 

12 4 y 



{qagGq) 



+ 



f rus mq\ 
V12 4 7 



{sagGs) ^ rnqTUs 



8q^ 



^G' 

TT 



16 



Snag 
2V 



(A.5) 



[m 



2 2\ 

Trig) 



727r3 



- + ln-^ 



+ 25 In — 



6g 



TT 



In- 



lOniqas , , 



6g4 



TT 



{sagGs) - -^{qagGq) 



1 , 



1 - - In 



+ 



in a. 



27g^ 



[{qqf-{ssf]. 



(A.6) 



The expression for U^^ has already been given in Ref. [5], together with Hq^^^ 2)' leading 
order in SU(3) breaking. The terms in and niq are new. The above expressions fulfill 
the relation (A. 4). 

At this point a few comments are in order concerning the structure of these formulas. 
The reader may have noticed that the Wilson coefficient of the gluon condensate contribu- 
tions to the above correlation functions contain infrared sensitive terms ~ ln(— g^/mg ^). 
These terms appear to violate the structure of the operator product expansion which stipu- 
lates that long- and short-distance contributions be properly factorised and all long-distance 
contributions be absorbed into the condensates, leaving purely short-distance Wilson coef- 
ficients which must be analytic in rriq^s- As discussed in Ref. [19], the appearance of terms 
logarithmic in ruq^s is due to the fact that the above expressions arc obtained using Wick's 
theorem to calculate the condensate contributions, which implies that the condensates arc 
normal-ordered: (O) = (0| : O : |0). Recasting the OPE in terms of non-normal-ordered 
operators, all infrared sensitive terms can be absorbed into the corresponding condensates. 
Indeed, using [19] 



a. 



{0\sgGs\0) = (0| : sgGs : |0) + log -^{0\:^G': |0), 



TT 



and the corresponding formula for q quarks, all terms ~ In ^ can be absorbed into the 
mixed quark-quark-gluon condensate and the resulting Wilson-coefficients can be expanded 
in powers of rrig g. In calculating the sum rules, we hence will use 



In 



In 



m 



q,s 
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As for the structure of the ultraviolet logarithms ~ ln(— they follow from the mixing 
of the gluonic operator q{gGafj,)ij^j5S with various quark-bilinear operators as given in 
Eq. (20) in Ref. [5]. 



A.2 Kl(K*) 

The correlation functions used to determine k\{K*) are very similar to those in the previous 
subsection. We choose = ij^ and T2 = CFp^ to obtain the pure-parity correlation function 

and r2 = 7/j which yields 

n^") and receive contributions from 1 states only and U-q\ from both 1 and O"*" 
states. Another possible choice is r2 = 1 which yields the pure-parity correlation function 

n„(g) = g.ng)(g2) (A.9) 

with contributions from only 0+ states. 11^'' and n^^2) related by the equation of 
motion for the vector current: 

ng^(g^) + q^u'a^^iQ^) = {tjIs - mg)u'a^ + contact terms. (A.IO) 

The expression for 11^^ was given in Ref. [5], the other correlation functions are obtained 
by the simple replacements 

■mq^-mq,{qq)^-{qq),{qagGq)^-{qagGq) 

ng)(g2) = ng)(g2) (A. 12) 

mq^-mq, {qq)^-{qq), {q(rgGq)^-(qagGq) 

which follows from the chiral structure of the correlation functions. 



A.3 K^iK*) 

For K,^{K*), r1 is given by a^^ and for r2 we choose a^s- The resulting correlation function 
has contributions from both 1" and 1+ states and can be written as 



where the projectors P^^ are given by 



[{ga-rQuQs - {a ^ /?}) - ({7 ^ <^})] , 



(A.13) 



3,a/37<5 



[{ga^qpqs + {a^ 13}) - ({7 ^ 5})] , 



q- 
1 



4,a/37<5 
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projects onto the twist-3 matrix element k^{K*), Pi onto K,j^{K*) and P^^ onto the 
contribution from 1"*" intermediate states. As usual, Hap-yS must not have a singularity at 
— which implies 

n&(o) + ng4(o) = o. 

That means that one can construct a mixed-parity sum rule from IIg,a = i^GAi^l'^) + 

4(0.'^)) / (f which has lower dimension than the pure-parity sum rule obtained from 
alone. We find^ 



^GAil') = (ml - ml) 
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-\- 2mqms In — 
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(A.14) 
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31n2^-171n^ 



+ 



+ 



Stt 
1 



.2 4 



(-In ^ + ;^ ) {qq) 



6 V /^^ 3 



37r 



12g2 



{qagGq){2ms - niq) 



I2q' 



{s(TgGs){ms — 2mq) 



.2 4 



(-In -1- + ;: ) (55) 



IJ? 3 



+ 



1 /"■■'■, 



24g2 \ TT 
+ 0-((g-g)2-(s.)2), 



2 - In 



2 - In- 



m 

- 2mqms In — ^ 



(A.15) 



nG,4(g') = {ml -ml) 



a. 



727r3 



31n2-^ - 141n 
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+ 



9V 



[6(ss) - (gg)] - [6(gg) 



97rg2 



+ 



1 

6^ 



[mq{qagGq) - ms{sagGs)) 



^We also give g^-independent terms in the quark condensate contribution to Hq'^ because they are 
needed for calculating IIg,4- Note that for IIq^ these terms are affected by finite counterterms as discussed 
in Ref. [6], which however cancel in the sum Hq ^ + U.^^^. 
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B Borel Transforms 



rrig 



+ ml 



2-ln^ 



(A.16) 



QCD sum rules are obtained from the Borel transforms of the correlation functions listed 
in the previous section. Most of the transforms are straightforward, except for those of 
expressions of type l/(g2)'^ln(— g^//^^), which can, however, be conveniently calculated 
using the formula 

TT i (— ctji (1 + a) 

with s — —(f. We then obtain, including continuum subtraction of contributions from 

s > So, 
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